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One key diagnostic discriminating among physical emis- 
sion mechanisms is the motion of the X-ray source: hydrody- 
namical numerical models of continuous proto stellar jets plow- 
ing through a uniform medium show an X-ray emitting shock 
front moving at several hundreds km/s. In the nearest X-ray 
emitting protostellar jet, HH 154, this is detectable, with the 
spatial resolution of the Chandra X-ray observatory, over a few 
years baseline, allowing a robust discrimination among differ- 
ent mechanisms. We have performed, in October 2005, a deep 
Chandra X-ray observation of HH 154. Comparison with the 
previous (2001) Chandra observation allows to detect proper 
motion down to the level predicted by models of X-ray emitting 
shocks in the jet. The 2005 Chandra observation of HH 154 
\^ shows unexpected morphological changes of the X-ray emis- 
sion in comparison with the 2001 data. Two components are 
. present: a stronger, point-like component with no detectable 

■ motion and a weaker component which has expanded in size by 
^ approximately 300 AU over the 4 years time base of the two ob- 

u servations. This expansion corresponds to approximately 500 
, km/s, very close to the velocity of the X-ray emitting shock in 

QQ the simple theoretical models. The 2005 data show a more com- 
, plex system than initially thought (and modeled), with multi- 

■ pie components with different properties. The observed mor- 
^ phology is possibly indicating a pulsed jet propagating through 

a non-homogeneous medium, likely with medium density de- 
creasing with distance from the driving source. Detailed theo- 
ff^ retical modeling and deeper X-ray observations will be needed 
, to understand the physics of this fascinating class of sources. 
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Abstract. Protostellar jets are a new class of X-ray sources which has been discovered with both XMM-Newton and Chandra. 
The mechanism responsible for the X-ray emission is still not clear. Self-shocking in jets, shocks where the jet hits the sur- 
rounding medium, reflected or scattered stellar X-ray emission have all been invoked as possible explanations. 



1. Introduction 

Soft X-ray emission from jets driven by protostellar sources 
has been discovered in 1999 thanks to the improved sensitiv- 
ity and spatial resolution oflTered by both the XMM-Newton 
and Chandra obs ervatories. HH2 in Orion has been the first 
source detected ( Prav do et all200ll) . with HH 154 (the jet orig- 
inating from the embedded binary protostar IRS 5 in the LI 55 1 
star-forming region) fol lowing; shortly thereafter (Favata et al. 
l2002t lBallv^tal1l2003h . While more sources have been de- 
tected later, HH 154 remains the best source to study thanks 
to both its being nearby (d ^ 150 pc) and to its being observ- 
able through a moderate absorbing column density (Ay ^ 7 
mag) while the parent star IRS 5 is hidden by ca. 150 mag of 
absorption. This allows to study the X-ray emission from the 
jet with good detail and without any glare from the parent star. 

X-ray emission from HH 154 was observed with XMM- 
Newton (Favata et al. 2002), allowing to determine the temper- 
ature of the source at about 4 MK. The limited spatial reso- 
lution of XMM-Newton did not allow to pinpoint the source 
of X-ray emission, and Favata et al. (2002) speculated that the 
source of X-ray emission might be located within the jet's 
working surface (i.e. the shock formed when the jet hits the 
circumstellar material), located in 2005 at some 15 arcsec 
from the parent sta r. The spatial resolution of Chandra allowed 
iBallv et all ilOO^ to reject this hypothesis, showing that the 
X-ray emission is located very close to the protostar itself and 
thus to the jet's base. 

A number of hypotheses have been advanced re lative to the 
physic al mechanism causing the X-ray emission. iBally et alJ 
(I2OO3I) speculated on the various possibilities, including re- 
flected or scattered stellar X-ray emission, internal shocks in 
the jet and shocks formed when the jet hits the cavity wall. 



Later, iBonito et"all (l2004l) performed detailed hydrodynamic 
numerical modeling of continuous jets propagating in a uni- 
form medium, showing that the X-ray emission observed from 
HH 154 can be explained in such a scenario as coming from 
the shock forming where the jet interacts with ambient matter, 
with reasonable assumptions about e.g. densities and velocities 
for both the jet and the medium. 

Other proposed mechanisms for the origin of the observed 
X-ray emission predict diflTerent properties, which can be used 
as discriminating diagnostics. For example, scattered or re- 
flected stellar X-ray emission should show the same temporal 
and spectral variability of the emission from the parent star, in 
particular the ubiquitous flaring obs erved in PMS stars. O n the 
other hand, the numerical models of 'Bonit o et alJ (|2004') show 
the X-ray emitting shocks to be moving at speeds of ca. 500 
km/s, with little temporal variability. 

Th e motion of the X-ray emitting shock predicted by 
iBonito et al. ( 2004) 's models is sufi&ciently large as to be de- 
tectable with the high spatial resolution of Chandra over a tem- 
poral baseline of a few years. With the main aim of constraining 
the physical mechanism driving the X-ray emission from the 
jet we have performed, in 2005, a deep Chandra observation of 
HH 154. Its detailed comparison with the 2001 observation al- 
lows a strong test of whether the X-ray source indeed displays 
proper motion, and thus a test of whether the "moving shock" 
model is indeed consistent with data. 

2. Data analysis and results 

The Chandra observation used here was performed on 27 Oct. 
2005, for a total duration of 97 ks. The observation was cen- 
tered at 04:31:32.9, +18:08:04.7 (FK5) and was performed 
with the ACIS-I camera. We processed the observation using 
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Fig. 1. The top panel shows the X-ray source associated with the HH 154 jet, in the 2005 (left) and 2001 (right) Chandra ob- 
servations. The small white circle in the 2005 observation shows the position of the protostar driving the jet (LI 551 IRS5, at 
04:31:34.2 -hi 8:08:05). The scale of the image is 0. 5 arcsec/pixel, nor th is up, east left. While a hint of an extended component 
is present already in the 2001 data (as discussed by Bally et in i2003h . the 2005 observation clearly shows two components to 
be present, a point-like source which has the same position as in 2001, and an expanding component which moved, in the jet's 
direction and away from the driving star, by approximately 2 arcsec between 2001 and 2005. The color scale is the same for both 
observations. For comparison the bottom panels show (in the same color scale) a background AGN at 2.5 arcmin from HH 154, 
demonstrating in both cases the sharpness of the Chandra PSF for point-like sources. 



the CIAO package, and filtered the data to retain only photons 
in the 0.3 to 4.0 keV band. To allow an uniform comparison, 
we have al s o repr ocessed the 2001 observation (described by 
iBallv et all (l2003h ) using the same approach. Individual X-ray 
sources have been detected using the standard wavelet-based 
algorithm of CIAO. A small positional shift is present between 
the two observations, as shown by both the stellar and the ex- 
tragalactic X-ray sources in the field. To allow detection of the 
expected small proper motion of the jet X-ray source, we have 
registered the two images by applying a shift determined by 
minimizing the coordinate diff'erences of the 7 brightest X-ray 
sources in the field which are also members of LI 551 (to min- 
imize the eff'ects of proper motion). The shift between the two 
observations is 2.43 ± 0.28 arcsec, a value within the absolute 
attitude reconstruction accuracy expected for Chandra. 

The first surprising result of our 2005 observations is that 
the X-ray source associated with the jet is clearly extended, 
as shown in Fig.[l] which shows both the 2001 and the 2005 
Chandra observations of HH 154. The position of the driving 
protost ar (L1551 IRS 5) is also plotted, from the VLA observa- 
tions of lRodri^uez etalJ (l2003h . 

Already for the 2001 observation iBallv etalJ (l2003h had 
suggested that the jet X-ray source might show a hint of be- 
ing extended; in our 2005 observation not only do we confirm 
that the source indeed has an extended component, but we also 
show that the morphology has changed, with the size of the 
extended component clearly increasing with respect to 2001. 
The bottom panels of Fig. [l] show the image of a background 



AGN located at 2.5 arcmin from the jet source. In both 2001 
and 2005 the source is clearly point-like, and with a diff'erent 
morphology than the jet, showing the latter to be resolved by 
Chandra. 

In the 2005 observation the source appears to have both 
an unresolved, point-like component and an extended compo- 
nent. While the unresolved component does not show a de- 
tectable motion with respect to 2001, the extended component 
has clearly become more elongated along the jet's direction, as 
expected from a shock moving within the jet and away from the 
parent star. The total count rate from both sources is compatible 
with their being constant, with no short-term variability visi- 
ble. Fig. 12 shows a horizontal cut through HH 154 in the 2005 
observation, showing the relative intensity of the stationary and 
moving sources. 

To confirm the extended nature of the source in the 2005 
observation, we performed an adaptive smoothing of the data 
with a Gaussian kernel, using the software tool csmooth avail- 
able in the CIAO software. The results, shown in Fig. (3] clearly 
show the source associated with HH 154 to be elongated in 
2005. In 2001 the same source was compatible with being 
point- like. At the same time, the AGN as shown in Fig. [l] has 
a symmetric, point-like shape in both smoothed images (not 
shown in Fig. [3). This again demonstrates the change in mor- 
phology of the X-ray source associated with the jet over the 4 
years. 

Between the 2001 and 2005 observations the source extent 
has increased by 2 arcsec, which, at the distance of L1551, cor- 
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respond to 4.5 x 10^^ cm, or 300 AU. This results in a projected 
velocity of ca. 33 km/s, v^hich, assu ming the inclination of the 
jet to be 45 deg (iLiseau et al.ll2005l) . implies a shock velocity 
of ca. 460 km/s, or about 100 AU/yr. 

We have also performed a simple spectral analysis of the 
source's emission. The integrated source has 63 photons, which 
we have fit with a simple absorbed thermal plasma (apec model 
in XSPEC), resulting in a temperature of approximately 0.5 
keV and a column density of 1.5 x 10^^ cm"^. This is essen- 
tially identical to the integrated X-ray spectral characteristics 
determined with the Sep. 2000 XMM-Newton observation, and 
compatible with the Chandra 2001 observation, showing lim- 
ited spectral variability in the source, even in the presence of 
significant morphological evolution. The number of X-ray pho- 
tons present in the moving component of the source is too small 
to attempt a separate spectral analysis. While the energy dis- 
tribution of the photons from the moving source is compatible 
with the spectrum of the non-moving component, there is a hint 
of more photons present at lower energy, which would be com- 
patible either with a lower absorbing column density (i.e. with 
the moving component coming out of the large absorption re- 
gion in which the jet is born) or with a lower temperature. 

The count rate from the integrated source in the 2005 ob- 
servation is 0.65 cts/ks, compatible within the Poisson errors 
with the count rate observed in the Chandra 2001 observation 
(0.77 cts/ks), showing limited (if any) changes in the source's 
luminosity (given the lack of observed spectral evolution). 
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Fig. 2. A horizontal cut through the 2005 Chandra HH 154 im- 
age, obtained by adding the number of photons in the two rows 
where the stationary source peaks, showing the relative inten- 
sity of the stationary and the moving components. 



Fig. 3. The smoothed image of the X-ray sources associated 
with HH 154 are shown, in the left panel for the 2005 data and 
in the right panel for the 2001 data. The 2005 data show the 
presence of an elongated component not present in the 2001 
data. 



3. Modeling implications 

The scenario of lBonito et al.l(l2004 l). which the observation dis- 
cussed here was designed to test, made clear predictions about 
the evolution of the X-ray emission due to the jet's shock. In 
particular, in the model which best reproduced HH 154's X- 
ray luminosity and spectrum the X-ray source was point-like, 
located at the head of the jet and had a continuous motion of 
about 500 km/s. 

The observational results discussed here are more complex 
than expected: the X-ray source displays both an apparently 
point-like component, with no measured motion, and a prob- 
ably extended component which is moving (from the diff'er- 
ence in position between the 2001 and 2005 observations) at 
about 500 km/s, a value very close to the value predicted by the 
Bonito et al. (2004) model. However, the apparently stationary 
compo nent is not in agreement with the simple iBonito et all 
(l2004 model. 

The assumptions behind this model were simple, aiming at 
elucidating the basic physics. In particular, the jet is assumed 
to be continuously powered, with constant density of inflow- 
ing material, and the ambient medium through which the jet 
propagates is assumed to be uniform. Neither of these two as- 
sumptions is likely to be true in realistic conditions. 

The optical observations of HH 154 in the last two decades 
show evidence of the jet being highly inhomogeneous, with 
density varying by some tw o orders of magnitude within the jet 
itself (Fridlun d et al.l2005l) . so that the jet is likely to be driven 
in an eruptive, rather than continuous, form, with new visible 
features appearing on a time scale of few years. Evidence for 
the jet to be pulsed, rather than continuous, comes also from the 
decreasing velocity of individual blobs through the jet: while 
the de-projected velocity from optical spectroscopy is up to 600 
km/s close to the star, at 3 arcsec from it it begins to drop, de- 
creasing down to 160 km/s for the working surface which is 
located at some 15 arcsec from the star in the 2005 observa- 
tions - see Fridlund et al. (2005) for details. In the presence of 
a continuous driving of the jet the shock front would be ex- 
pected to travel at constant velocity all the way to 15 arcsec; 
the observed velocity decrease is thus likely indicating a sig- 
nificant pulsating component. 

At the same time, the assumption of a homogeneous 
medium through which the jet propagates is unlikely to be 
correct: the high level of structuring of the medium is shown 
for example by the fastly varying circumstellar absorption (the 
driving star being hidden by some 150 mag of material, while 
the jet has less than 10 mag of absorption). Even along the 
jet optical observations show the absorption varying between 
2 and at least 6 mag (Fridlund et al. 2005). 

Although the predictions of the simple model of a continu- 
ous jet ramming into a homogeneous medium of 'Bonit o et alJ 
12004) is not compatible with these observational results, the 
high level of structuring in both the jet and the medium can 
possibly explain the phenomenology observed in X-rays, with 
the same basic physics on which the simple models are based. 
The density of circumstellar material is very likely to decrease 
significantly with distance from the jet-driving source (as ob- 
served in the structure of the absorbing material), and thus the 
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individual jet pulses (or "bullets") will initially meet a dense 
me dium, and then later propagate through a thinner medium. 

iBonito et"al] ('2006) have carri ed ou t, building on the ap- 
proach first used by Bonitoetal] (12004 ). an extensive explo- 
ration of the parameter space of both density contrast and ve- 
locity for a uniform, constantly driven jet. One result is that the 
X-ray luminosity of the shock front depends on both the density 
contrast and the Mach number with which the jet plows through 
the medium. In the presence of a decreasing ambient density 
and increasing temperature (due to pressure equilibrium), the 
jet's luminosity will decrease as it plows though the medium. 
In addition, in the case of an isolated traveling plasma bullet, 
the temperature of the X-ray emitting plasma is expected to 
gradually decrease due to radiative and conductive losses. 

One possible scenario explaining the observed phe- 
nomenology is therefore that the driving source is continuously 
injecting fresh blobs of material into the jet. As each blob rams 
into the high density medium near the star it will cause the 
X-ray emission seen as the stationary X-ray source closer to 
the protostar. As a given blob plows through the high density 
medium and starts to move through the thinner medium farther 
away from the star, its X-ray luminosity will decrease while 
it moves away, forming the fainter moving X-ray source well 
evident in our 2005 Chandra observation. 

This scenario is of course schematic, and it will require 
detailed modeling before it can be tested in detail. I t is how- 
ever su££orted_b^_optical dpridlund & Liseaulll994l) and ra- 
dio dRodrfguez et al.lEo03l) observations of new blobs appear- 
ing repeatedly from behind the obscured region every few 
years. The basic characteristics of our scenario are based on 
the detailed s tudy o f the simple, homogeneous jets studied by 
IBonito et"al] (l2006l) and rely on the same physics. Proper test- 
ing of the scenario (a pulsed jet in a non-homogeneous circum- 
stellar medium) will also require more sensitive X-ray observa- 
tions, with sufficient statistics as to allow a determination of the 
spectrum of the moving component (the one plowing through 
the thinning circumstellar medium), whose temperature is ex- 
pected to decrease with distance from the star. 

4. Discussion 

The Chandra 2005 observation discussed here was designed 
to provide a clear test of whether the X-ray emission from 
HH 154 is due to moving shocks rather than stationary mech- 
anisms such as scattered stellar X-rays. The motion of the X- 
ray source unambiguously shows that moving shocks play a 
significant role in the process. However, the observed X-ray 
morphology is more complex than what would be observed in 
the presence of a continuously driven jet propagating through 
a uniform medium. This is not surprising given the evidence 
from optical observation of a pulsating jet driven from an erup- 
tive process, and the observed non homogeneous distribution 
of absorbing material in the region surrounding the jet. 

Detailed study and modeling of a number of scenarios 
along the above lines will be needed to determine in de- 
tail which conditions can reproduce the observed morphology. 
While this work is ongoing, simple considerations based on our 
current models of uniform jets lead us to propose that the most 



likely scenario explaining the observations is one of an eruptive 
jet propagating within a medium with strong density gradients. 
In addition to the (already ongoing) detailed modeling, new, 
more sensitive X-ray observations will be needed in the future, 
for example to determine the temperature of the plasma associ- 
ated with the moving X-ray source, to ascertain whether this is, 
as expected in such a model, lower than the X-ray temperature 
associated with the stationary component. These measurement 
would also provide a direct determination of the scale of den- 
sity changes in the circumstellar medium around the star, pro- 
viding an unique probe of such material. Also, they will provide 
an unique probe of the high energy phenomena at the base of 
the jet, close to the launch site, constraining the physics of the 
jet formation. 
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